Morphology of two dimensional fracture surfaces 
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Abstract. We consider the morphology of two dimensional cracks observed in experimental 
results obtained from paper samples and compare these results with the numerical simulations 
of the random fuse model (RFM). We demonstrate that the data obey multiscaling at small 
scales but cross over to self-affine scaling at larger scales. Next, we show that the roughness 
exponent of the random fuse model is recovered by a simpler model that produces a connected 
crack, while a directed crack yields a different result, close to a random walk. We discuss the 
multiscaling behavior of all these models. 
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Understanding the statistical properties of fracture surfaces has been an important 
theoretical challenge for the past twenty years, starting from the first pioneering experimental 
evidence of self-affinity provided by Mandelbrot et al. (1). Most experimental results reported 
in the past for three dimensional fracture surfaces suggested the presence of a universal 
roughness exponent in the range £ ~ 0.75 — 0.85 (0). The scaling regime is sometimes 
quite impressive, spanning five decades in metallic alloys (|2J). Here, it was argued that a 
different exponent (i.e. £ ~ 0.4 — 0.5) should describe the small scales, with a crossover 
originally interpreted as a dynamic effect. This exponent would correspond to the quasistatic 
limit while the large scale exponent to the effect of finite velocities (2). It was recently pointed 
out that the short-scale value is not present in silica glass, even when cracks move at extremely 
low velocities In addition, in granite and sandstone, one only observes the "small scale" 
exponent even at high velocities (4). The current interpretation associates the value C — 0.75 
with rupture processes occurring inside the fracture process zone, where elastic interactions 
would be screened and the value ( ~ 0.45 with large scale elastic fracture <E). The authors 
of Refs. (0;E|) were also able to show that the fracture surface is anisotropic, with different 
exponents in parallel and perpendicular directions to the crack propagation. In addition, we 
have to remark that the measured roughness exponent describes only the local properties of the 
surface. The fracture surface in many cases exhibits anomalous scaling: the global exponent 
describing the scaling of the crack width with the sample size is larger than the local exponent 
measured on a single sample (refs. (0; Q; 0)). It is thus necessary to define two roughness 
exponents a global (Q and a local Q oc . 

Two dimensional fracture surfaces are in principle simpler to analyze than the three- 
dimensional surfaces which can be anisotropic and since the crack surface reduces to a line 
in 2D. The existing experimental results, obtained mainly in paper, point towards a (local) 
roughness exponent in the range £ ~ 0.6 — 0.7 (SIlO; 11; 12). However, even for ordinary, 
industrial paper itself there are numerous values available that are significantly higher than 
( = 0.7 (examples are found in Ref. (|13f) V It is not known at this time whether this variation 
in ( values is a reflection of difficulty in experimentally measuring £, or that it is not really 
universal but depends on material parameters such as ductility and anisotropy. Recently, 
Bouchbinder et al. have indicated a scenario in which the the crack line h(x) has more 
complicated structure, exhibiting multiscaling behavior. This implies a non-constant exponent 
ot q , for the g-th order correlation function C q (x) = (\h(x + y) — h(y)\ q ) 1 ' q ~ x aq , Q. This 
result would strongly put into question the existence of a well defined roughness exponent in 
two dimensional fracture. It should also be noted that so far there is no experimental evidence 
for the presence of anomalous scaling in two dimensions. 

From the theoretical point of view, two dimensional fracture could appear as a relatively 
simple problem if we consider the crack surface as the trace left by a point (i.e. the crack 
tip) moving through a disordered elastic medium ([l5; 16). A similar idealization is also used 
in three dimensions where the crack tip is replaced by a deforming crack line front. Under 
mode I quasistatic loading, the fracture surface is shown to be only logarithmically rough 
in three dimensions, in contrast with the experiments (|17|) . This suggests that to understand 
the experiments one should consider additional ingredients, such as damage nucleation ahead 
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of the tip ilat) . crack branching or elastodynamics effects. In this perspective, disordered 



jrspe 
,19; 



lattice models provide an alternative way to describe the phenomenon (19; 20). In these the 
elastic medium is described by a network of springs with random failure thresholds. In the 
simplest approximation of a scalar displacement, one recovers the random fuse model (RFM) 
where a lattice of fuses with random threshold are subject to an increasing external voltage 
(EH I22L The model has been numerically simulated to obtain the roughness of the fracture 
surface in two (|23; 24; 25; 26) and three dimensions J27; 24; 28; 29). In two dimensions, the 
roughness exponent Q oc ~ 0.7 is reasonably close to the experimental results. In addition, 
recent simulations reveal that anomalous scaling is also present, although the effect is small 
(i.e. C-Ctoc = 0.1) Q. 

In this letter we first point out that the multiscaling behavior observed in Ref. dl4l) is in 
fact a small scale effect due to the fluctuations induced by the fibrous structure of paper. To 
this end we re-analyze cracks produced in 6500 mm long paper samples (UJJ). We next turn 
our attention to the RFM and observe that similar corrections exist there as well. In order 
to understand the mechanism underlying the roughness in the RFM we introduce two simple 
models that describe the growth of a single connected crack. The numerical results show that 
a fracture process zone (FPZ) is essential in recovering the RFM scaling. We thus conclude 
that considering the motion of a crack tip represents an oversimplification of the problem, 
leading to quantitatively different results. 

Paper would seem to be a good test material for fracture surface analysis since it has 
a disordered structure and can be considered effectively two-dimensinal (for an in-depth 
discussion see Ref. (30)). However, as in any similar problem, the coarse-grained behavior 
can only be seen at scales that are larger than those associated with the microscopic details. 
Fig. d illustrates (left panel) the intricacies of fracture line analysis in paper. First of all, on 
a scale below (typically) 0. 1 mm, paper is in fact three-dimensional and the fracture surface 
is no longer a line. Second, the damage is diffusive on the typical FPZ scale, which ranges in 
industrial papers up to 2-3 mm. The greyscale in the figure shows in fact such a damage profile 
and illustrates that on such scales the final single- valued fracture line will have steep gradients 
Ah. On a scale Ax of the order of the micro structural details the distribution P(Ah)(Ax) 
displays non-Gaussian tails O). The C q (x) can be sensitive to such tails, when the height 
differences \h(x + x') — h(x')\ that dominate C q are also important for the q moment of 
P(Ah). This is demonstrated in the right-hand panel of Fig. d where we plot the C q for 
various moments q ranging from q = 1/4 to 6. It is clear that there is apparent multiscaling 
related to the P(Ah) (see (UJJ)) and that beyond a cross-over scale of a few mm one obtains 
self-affine scaling with ( = 0.64 in agreement with other measures (fill) . In passing, we note 
that the structure of paper exhibits correlations (due to so-called floes) over still larger scales 
than the FPZ size extending up to several mm. The cross-over scale might be related to these 
correlations. 

In the RFM (|2l|), one considers a triangular lattice with fuses having all the same 
conductance and random breaking thresholds t, uniformly distributed between and 1 . The 
burning of a fuse occurs irreversibly, whenever the electrical current in the fuse exceeds 
its threshold t. Periodic boundary conditions are imposed in the horizontal direction and a 
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Figure 1. The diffuse crack surface in paper and the thresholded "interface". At small 
scale it is not obvious how to define a single valued interface and jumps associated with 
overhangs are unavoidable (left). As a result of this one sees an apparent multiscaling at small 
scales, followed by self affine scaling at larger scales. The curves represent the q— correlation 
functions, normalized by their maximum value for clarity (right). The dashed line has an 
exponent £ = 0.64. 

constant voltage difference, V, is applied between the top and the bottom of lattice system 
bus bars. Numerically, a unit voltage difference, V = 1, is set between the bus bars and 
the Kirchhoff equations are solved to determine the currents in the fuses. Subsequently, for 
each fuse j, the ratio between the current ij and the breaking threshold tj is evaluated, and 
the fuse j c having the largest value, max^, is irreversibly removed (burnt). The current is 
redistributed instantaneously after a fuse is burnt implying that the current relaxation in the 
lattice system is much faster than the breaking of a fuse. Each time a fuse is burnt, it is 
necessary to re-calculate the current redistribution in the lattice to determine the subsequent 
breaking of a fuse. The process of breaking of a fuse, one at a time, is repeated until the lattice 
system fails completely, producing an irregular fracture surface. 

Here, we propose two variations to the RFM. (i) In the first variation, we impose that 
failure events form a connected crack, excluding damage nucleation in the bulk. This means 
that after breaking the weakest fuse, successive failure events are only allowed on fuses that 
are connected to the crack. Otherwise, the rules of this simplified model strictly follow those 
of the usual RFM. In effect this rule implies that we have a FPZ which is constrained to a 
distance r = 1 from an evolving crack, (ii) In the second variation, in addition to the variation 
(i), we also do not allow for any crack branching, nor turning backwards. We thus only break 
one of the three fuses connected to the crack tips (i.e. the one with the largest current/threshold 
ratio). This leads to a single directed crack, with the surface being the trail left by the crack 
tip. 

We simulate the variations (i) and (ii) on triangular lattices of linear sizes L = 
128, 192, 256, 320, 512 with uniformly distributed disorder. The final crack in case (i) and 
in the RFM typically displays some limited amount of dangling ends and overhangs. These 
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are removed to obtain a single valued crack line h(x). This is not necessary for case (ii) since 
the crack is by definition single valued in this scenario. Several methods have been devised 
to characterize the roughness of a self-affine interface. The power spectrum is believed to 
provide one of the most reliable estimates of the roughness exponent, and decays as 

S(k) = J2e i{27Txk/L) (h(x)h(0)) ~ AT (2C+1) . (1) 



It was shown in Ref. (1261) that the RFM displays anomalous scaling, and in this case, the 
power spectrum scales as S(k, L) ~ fc-( 2 Cioc+i)£2(c-c !oc )_ 



In Fig. |3 we compare the power spectra of the cracks obtained from the RFM (1261) and 
from variations (i) and (ii). The RFM and variation (i) follow essentially the same scaling 
behavior with <^ oc ~ 0.7 and ( ~ 0.8. On the other hand, the directed crack scenario (variation 
(ii)) differs considerably: there is no anomalous scaling and the exponent is significantly lower 
(i.e. C — Ooc — 0.46). This result is very close to the random walk exponent ( = 1/2 that 
would be found in the directed crack model wherein the effects of current are completely 
ignored. In this scenario, the crack tip would move up and down depending only on the 
smaller threshold. The fact that the exponent is slightly smaller than £ = 1/2 could be 
a numerical artifact associated with random walk on triangular lattice topology. Hence, to 
recover the RFM roughness exponents, the presence of crack branching, as in variation (i), is 
essential. This result indicates that the presence of distant damage nucleation and coalescence 
with other cracks is irr 




Figure 2. The power spectrum of the crack S(k, L) obtained from different models for 
different lattice sizes in log-log scale. The slope defines the local exponent as —(2Q oc + 1). 
The spectra for all of the different lattice sizes can be collapsed using simple self-affine scaling 
only for directed cracks, yielding £ = Q oc = 0.46. Both the RFM and the single crack model 
show instead anomalous scaling with £; oc = 0.7 and £ = 0.8 



In Fig. 3, we analyze the multscaling behavior of the RFM and its two variations. As 
expected, the directed crack does not display any multiscaling since this is akin to a random 
walk. On the other hand, multiscaling is observed both in RFM and in variation (i) models. 
We notice that deviations from simple self-affine scaling are stronger for high q values, while 
for low q one recovers the a q = Q oc . The origin of this behavior is related to the removal of 
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overhangs on the crack surface, a process that inevitably produces steps in the single valued 
crack profile. As discussed by Mitchell, adding random steps to a self-affine profile yields 
an apparent multiscaling over small scales and for q > 1, while for low q values, one obtains 
a q — Cioc (1321) . Therefore, we conclude that the apparent multiscaling in RFM and in variation 
(i) models at small scales and for q > 1 is due to the process of removal of overhangs, and that 
pure self-affinity is recovered at larger scales. This can be further elaborated by considering 
the height difference Ah(l) = \h(x+l)—h(x) \ . We have checked its scaling with / and varying 
L, and three conclusions can be drawn: (1) (Ah(l)) ~ fa oc as expected. (2) The L-scaling 
exhibits similar anomalous scaling as the width and power spectrum data. (3) The crack 
profile is stationary (e.g. (A/t(l)} does not depend on x). Hence, the observed anomalous 
scaling is not due to non-stationarity of the crack growth process but rather to an intrinsic 
dependence of the era 




Figure 3. The g-correlation functions, normalized by their maximum values, for directed 
cracks (top), connected cracks and RFM (bottom). Multiscaling is essentially absent for 
directed crack, while both connected craks and the RFM display deviations from simple 
scaling, especially at high q. Notice that the multiscaling behavior of connected cracks and 
RFM follows a very similar pattern. 



Finally Fig. 4 depicts the scaling of the histograms of Ah(l) for L = 512, and for 



REFERENCES 



7 



variation (i). Empirical analysis of 2d fracture data implies that for values of I such that 
the profiles are self-affine, the distributions of Ah often approach Gaussian ( dill) , see also 
for a similar, claim). Here, we can discern clearly a Gaussian central part, but the tails do 
not follow a Gaussian even for / = 64 (similar to what is seen in Ref. dill) ). Varying L reveals 
again the presence o 
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Figure 4. The (logarithms) of histograms of Ah(l) for I = 1, 8, 16, 32, 64 and L = 512. As a 
guide for the eye we present a Gaussian fit for I — 32. It can be seen that the central parts of 
the distributions are Gaussian, but for large Ah deviations exist, probably originating from the 
overhangs that are removed. In addition, we report the data for L = 256 and / = 16 to show 
the systematic deviations due to the anomalous scaling with L 



In summary, having analyzed the crack morphology of two dimensional fracture surfaces 
in experiments and models, we conclude that multiscaling is an artifact due to the removal of 
small scale overhangs and that self-affinity is recovered at large scales. The study of single 
crack variations of the original RFM indicates that the roughness properties of the fracture 
surface are due to the damage accumulation within the FPZ surrounding the crack, whereas 
the diffusive damage nucleation distributed homogeneously over the rest of the geometry is 
irrelevant. Finally, a model of a moving crack tip appears to be an oversimplifaction of the 
problem, and yields quantitatively different results from those of the RFM. It is intriguing 
to speculate about the three dimensional case, where crack front line models have enjoyed 
a wide appeal in the literature but are always in quantitative disagreement with experiments. 
Simulations of connected cracks with large scale three dimensional fuse models could help 
clarify this issue. 
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